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1 Introduction 

We consider the initial boundary- value problem for the incompressible Navier- 
Stokes equations in half-space = {x^ > 0}: 

dWv = / mM^x(0,oo), (1.1) 
vi-,t)Ui = fort>0, (1.2) 

v{-,0) = vo in R\ (1.3) 

where vq is a sufficiently regular div-free field in M.^ with sufficiently fast 
decay for x — )■ oo and folaiRS = . 

Our main goal is to understand connections between possible blow-up of 
strong solutions and Liouville theorem for bounded ancient mild solutions, 
in the spirit considered for all space in [5]. We recall that the local- in-time 
existence of strong solutions was proved in all space in [S] and for bounded 
domains similar results appeared in [4]. As to unbounded domains, we refer, 
for example, to [7] or [2j. 

Let us consider the local-in-time strong solution v on its maximal interval 
of existence [0,T), which will be assumed to be finiteQ The time T is then 
the blow-up time. Such solutions are known to be unique (unlike the weak 
solutions, which on the other hand are global). We set 

g{t) = sup M(r) , (1.4) 

0<T<t 



"'^At the time of this writing it is unknown whether such solutions exist. 
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where 



M{t) = sup \v{x,t)\ . 



It is well-known that g(t) — t- +oo as t — )■ T_, see for example [S]. 

Our considerations are motivated by methods used in the theory of ge- 
ometric flows and based on scale-invariant transformations of the solution 
V when time is approaching T, see [5] for a more detailed discussion and 
references. In the case of the Navier-Stokes equations, the scale-invariant 
transformations have the form v{x, t), q{x, t) — )■ Af (Ax, XH) , A^g(Ax, XH). In 
the whole space, taking limits of sequences of suitably scaled solutions pro- 
duces so-called bounded ancient (backward) solutions to the Navier-Stokes 
equations. They are defined on the semi-infinite time interval ] — oo, 0] (back- 
ward in time), are bounded in L°° and in fact belong to a subclass called mild 
bounded ancient solutions. In [3], it has been shown that mild bounded an- 
cient solutions are infinitely smooth both in space and in time. Moreover, 
given that T is a blowup time, the limiting mild bounded ancient solution 
cannot vanish. It has been conjectured in [5] that any mild bounded ancient 
solution is a constant. This would rule out blow-ups of Type I in the case 
of the Cauchy problem (whole space). Let us recall that Type I blowup is 
usually defined by the inequality 



for any x G M'^ and any t < T, although more general definitions are also 
possible. 

To show a difference between bounded ancient solutions in the whole 
space and in the half space, we recall their definitions. We say that m is a 
bounded ancient solution of the Navier-Stokes equations if u is bounded in 
Q_ = M^x] — oo, 0[ and satisfies these equations in the sense of distributions 
with divergence free test functions, i.e.. 



v{x,t)\ < 



c 



(1.5) 



J yu ■ {dt(p + Aif) + u®u : V^pjdz = 
Q- 



(1.6) 



for any ^ e C^oiQ.) := G C^{Q-) : divy. = 0}; 




(1.7) 



Q- 
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for any q G C^{Q^). By scaling, we may assume that \u\ is bounded by one. 
There are simple nontrivial bounded ancient solutions of the form 

u{x,t) = a{t) (1.8) 

where a is an arbitrary bounded function of t only. 

A vector field u is called a mild bounded ancient solution if u is a bounded 
ancient solution and there exists a pressure field p e Loo (—00, 0; BMOiM.^)) 
such that 

J {u ■ {dt'^ + Aif) + u® u -.Vip^dz = — j pAiMifdz (1.9) 



for any G C^{Q-). It is not so difficult to see that any solution of the 
form (11.81) is a mild bounded ancient solution if and only if a{t) = constant. 
As already mentioned above, mild bounded ancient solutions are infinitely 
smooth and have the following property: for any A < 0, they can be presented 
in the form 

Ui{x,t) = T{x-y,t- A)ui{y,A)dy+ 




+ Kijm{x,y,t - T)uj{y,T)um{y,T)dydT, (1.10) 



where F is the well known heat kernel and K is obtained from the Oseen 
tensor by differentiation in spatial variables, see details in [13] and [5]. That 
is why those solutions are called mild ones. By the way, it is the original 
definition of mild bounded ancient solutions given in [5] and the equivalent 
definition in terms of the pressure appeared in [11] later on. 

The case of the half space is more complicated. First, it is not imme- 
diately clear how to understand the homogeneous Dirichlet boundary con- 
ditions when the velocity u is only bounded in Qt := {z = (x, t) : x G 
]R'j_, —00 < t < 0}. We shall use the following weak definition of (11. 2p : 

J (u- {dt^ + A(p) +u'S)u:V(p^dz = (1.11) 
Qt 
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for any ip G C^q{QS) with (/)(x',0,t) = for any x' G and for any 
-oo < t < 0; 

Ju-Vqdz = (1.12) 

for any g G C^((5-). 

In [3j, a class of simple non-trivial bounded ancient solutions to the 
Navier-Stokes equations has been presented. They describe a bounded shear 
flow in the half space and have the form 

U{x,t) = {ui{x3,t),U2{x3,t),0). (1.13) 

It has also been shown in [3] that there are no other non-trivial solutions to 
the linear Stokes system in the half space. It is unknown whether or not this 
is true for the nonlinear case. 

We now define mild bounded ancient solutions in a half space. 

Definition 1.1. A bounded function u is a mild bounded ancient solution if 
and only if there exists a pressure p such that p = p^ + p^ , where the even 
extension of p^ to the whole with respect to x^ is Loo{—oo, 0; BMO{M.^))- 
function, 

Ap^ = -divdivu (g) u (1-14) 

in with p\{x',0,t) = and p'^{-,t) is a harmonic function in M.^ whose 
gradient satisfies the estimate 

|Vp2(x,t)| < cln(2 + I/X3) (1.15) 

for all {x, t) G and has the property 

sup \Vp^{x,t)\ (1.16) 

as X3 ^ 00; u and p satisfy U.l^) and 

u ■ (dt^p + Ay?) + U0U : Vp + pdivLp^dxdt = (1-17) 

for any Lp G C^{Q-) with p{x', 0,t) = for any x' G and for any t < 0. 
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Remark 1.2. If u is a mild bounded ancient solution, then Vu G Loo{Q~!l). 
The function u is infinitely smooth in spatial variables in upper half space 
X3 > 0. 

As in the case of the whole space, mild bounded ancient solutions can 
be defined by formula fll.lOp . in which is replaced with Rj^ and kernels 
with their half-space analogues, for details see Section 2. The corresponding 
statement might be called the equivalence theorem. The proof of such a 
result is more involved than its whole space version and will be published 
elsewhere. 

It is worth noticing that a nontrivial solution of the form fll.lSp is a mild 
bounded ancient solutions if and only m = 0. Indeed, for solutions fll.l3p . the 
gradient of the pressure is a function of time only. By the above theorem, 
this is only possible if the gradient of the pressure is equal to zero. In turn, 
this means that each component u^, a = 1, 2 is a bounded ancient solution 
to the heat equation in {x^ > 0}x] — oo,0[ with the boundary condition 
Ua{0,t) = 0, which implies that Ua = 0. 

We believe that the following is true: 

Conjecture 1.3. There is no non-trivial mild bounded ancient solution to 
the Navier-Stokes equations in the half space. 

The validity of Conjecture 11.31 and the conjecture, made in [5] and men- 
tioned above, would rule out Type I blowups in the broad sense, when under- 
stood as blow-up solutions with a suitable bounded scale-invariant quantity. 

Let us state our main result. 

Theorem 1.4. Assume that initial boundary value problem U.l\) - in~^) has a 
solution that blows up at time T. There exists at least one non-trivial (non- 
zero) mild bounded ancient solution either in the whole space or in the half 
space. 

The appearance of mild bounded ancient solutions in the whole space is 
not surprising and it should be expected if solution v to original problem 
f ll.ip - fll.3p is smooth near the boundary x^ = near blowup time T. This 
scenario of the blowup can be called interior blowup. All other mild bounded 
ancient solutions then are related to the boundary blowup. We note that the 
boundary blowup could still lead to a mild bounded ancient solution in the 
whole space. This might happen when the velocity tends to infinity fast in 
comparison with the rate at which x is approaching the boundary. 
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An interesting consequence of the above theorem is the following state- 
ment. 

Proposition 1.5. There exists £ > (independent ofv) such that if 

\v{x,t)\<— (1.18) 

for all X e and t e]0, T[, the solution v does not blow up. 
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2 Prelimaries 



Given A < 0, consider the following Stokes problem in half space 
dtu — Au + Vp = — div F, div u = 

in x]A,0[, 

M(x',0,t) = 

for {x',t) e M2x]A,0[, 

u{x, A) = uo{x) 

for X e M^. 

In addition, we assume that uq is a divergence free and F and its first 
derivatives vanish at the boundary x^ = 0. If F and Uq are sufficiently smooth 
and decay sufficiently fast at infinity, a solution of the above problem can be 
presented in the following way, see 



u = u^ + u'^, (2.1) 

where 

u\x, t) = j G{x, y, t - A)uo{y)dy (2.2) 

for any x G and for any t g]A, 0[ and 

t 

u^(x,t) = - G{x,y,t-T)dwH{y,T)dydT (2.3) 



A 



for the same x and t, where H = F + p^I and p^ is a solution of the following 
Neumann boundary value problem 

Ap^ = — divdivF 

in and p\ = at 0:3 = 0. Here, G is the Green function for the Stokes 
system in the half space that has been studied in and has the form 

G = G^ + G\ (2.4) 

where 

Glj{x,y,t) = 6ij(T{x-y,t) - T{x-y*,t)^, 



X3 

d f f dE 

G%{x,y,t) = 4-^ —{x-z)r{z-y*,t)dz, G^3(x, t) = 0. 

IR2 



As in the case of the whole space, it is convenient to present the function 
in the following equivalent way 



t 

,2 



uf{x,t) = j J Kijm{x,y,t - T)Fjm{y,r)dydT (2.5) 



A ir3 

for the same x and t as above. The kernel K has been introduced in [13] and 
has the following structure 

■^ismix^ t) K igYTiix ^ t) -|- KismiXi Z^ t), (2.6) 

where Kism{x,z,t) is a linear combination of the terms 

dGij 

-7r^{x,z,t) 
dzk 

and Kismix, z, t) is a linear combination of the terms 

- — — / Gij{x,y,t)— {y,z)dy. 

dXaOXfS J oys 

Here, N^'^\x,y) = E{x — y) ± E{x — y*) with y* = {y', —y^) and E{x) is 
fundamental solution to the Laplace equation in R^. 

Let us outline how the above transformations can be done. Our arguments 
slightly differ from those that have been used in mentioned papers [B] and 
[13]. Consider the following boundary value problems 

Ay^rnnix, y, t) = G„„(x, y, t) (2.7) 

with d^mn/dy3{x,y,t) = if n < 3 and with ^jnn{x,y,t) = if n = 3 at 
1/3 = 0. And then integration by parts gives 
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The above splitting for the potential K can be easily derived with help of 
the initial boundary value problems for function see (12 .Tp . and some ele- 
mentary properties of Green functions G, A^~^, and N~ . 

The following estimates for G* and K have been obtained in papers [T2] . 
[13], and 



Q\a\+\l\Q2 



{x,y,t-A) < c{a,-f){t - A)-^ {t - A + xly^ X 



x{\x-y Y \ 2 exp 

where a' = (cti, a2), 7' = (71, 72), and I7I = 0, 1, 



t- A 



dG] 



(|x-l/|2 + t)2- 



d\G\x,y,t) 



< 



exp 



(2.8) 

(2.9) 
(2.10) 



t^{\x' -y'\'^ + xl + y^ + t)2 
for / = 0, 1. 

Let and be generated by and respectively. In particular, 
we need the following estimate 



\K'{x,y,t)\ < 



{\x - y*\^ + ty 



(2.11) 



which can be obtained by an elementary modification of arguments used in 
the proof of Proposition 3.1 in [13]. 
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3 Scaling 

It is not difficult to show that there exists a sequence {x^^\tk) with x^^^ E 
such that tk ^ T — and 

g{tk) = M{tk) = |t;(xW,tfc)| ^ oo. (3.1) 

Remark 3.1. From /i/ and from J^, it follows that sequence x'^^^ is bounded. 

We let Mk = M{tk). There are two main scenarios. In the ffist of the 
them, 

xi^^Mk oo (3.2) 

and we scale v and q so that 

u^''\y, s) = -^v{x, t), pk{y, s) = ig(a;, t), (3.3) 

where 

y = Mk{x-x^''^), s = Ml{t-tk). (3.4) 
By the above scaling, (II. ip and (11.21) are transformed into 

+ M^'^) ■ Vm(^) - AuW = -Vpfc, divM('=) = (3.5) 

in := {y = {y'.y,) : y' G R^, 1/3 > -xf M^jx] - 4M|,0[, 

u^^\y\-xfMk^) = Q (3.6) 

for any y' G and s g] — tfeM|,0[. And, moreover, according to (13. ip . we 
have 

|m('=)(0)| = 1. (3.7) 

In the second scenario, 

xf^Mfc ^ a G [0,oo[. (3.8) 

This suggests the same scaling (13.31) but with slightly different change of 
variables 

y' = Mk{x -x'^^^), y3 = MkX3, s = Ml{t-tk), (3.9) 
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In this case, fll.ip and fll.2p are transformed into system fl3.5p . which is vahd 
in R'^x] — tkM^, 0[, and into the boundary condition 

u^'\y',0,t)=0 (3.10) 

for any y' G and s g] — tkM^, 0[. Condition (13. 7p is replaced with 

|u('=)(0,4')m,,0)| = 1. (3.11) 

Our aim is to understand what happens if — )■ cxo. 

Without loss of generality, we may assume that the following statements 
are true: 

Scenario 1 There exists a divergence free function u G Loo(Q-) such that 
\u\ < 1 a.e. in Q_ and, for any R> 0, 

J u^^^ ■ wdxdt y" ■ wdxdt (3.12) 

Q{R) Q{R) 

for any w G Li{Q{R)). Here, Q{R) = B{R)x] - R\0[; 

Scenario 2 There exists a divergence free function u G Loo(Qt) such that 

|m| < 1 a.e. in and, for any R> 0, 

J u'^'^^ ■ wdxdt j u- wdxdt (3.13) 

Q+{R) Q+{R) 

for any w G Li{Q+{R)). Here, Q+{R) = B+{R)x] - R^,0[. 
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4 Scenario 1 

Our goal is to show that the hmit function u from (13.121) must be a mild 
bounded ancient solution in the whole space and satisfy the condition 

|m(0,0)| = 1. 

Fix an arbitrary A < 0. Let k is sufficiently large so that A > ~tkM^. 
We split w := u^'^^ into two parts 

w = + w'^ 

so that 

dtw^ - Aw^ + Vr^ = 0, div = 
in = X {?/3 > ~dk}y<]A, 0[, where dk = xf^ Mk oo, 

w^{y\ -dk,t) = 

for {y',t) eR'^x]A,0[, and 

w\y,A) = w{y,A) 

for ?/ G X {?/3 > -4}. 

The second part of u; is a solution to the following initial boundary value 
problem 

dtw'^ - Aw^ + Vr^ = div F + Wp^^^^ , div w"^ = Q 

in Q\, 

w'^{y\ -dk,t) = 

for {y',t) eR^x]A,0[, and 

w'{y,A) = 

for ?/ G X {?/3 > —dk}. Here, F = w ^ w and p^'-'^-* is defined by the 
following Newmann boundary value problem 

Ap'^'\y,t) = -divdivF(i/,t), {y,t) G Q^, 

pf\y',-dk,t) = {y',t)eR^x]A,0[. 
Keeping in mind (12 .4^ . we can present in the form 
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where 



w^''^{y,t)= j G^{y + dkCs^z + dke3,t - A)w{z,A)dz 



K2x{j/3>-dfc} 

and 



w^''^{y,t)= j G'^{y + dke3,z + dke3,t-A)w{z,A)dz. 

M.^x{y3>~dk} 

Elementary calculations, estimate fl2.8p . and the fact \w{y, A)\ < 1 ensure an 
upper bound for w^'^ 

|^^.^(y,t)|<c^^. (4.1) 
1/3 + dk 

The similar decomposition can be exploited in order to evaluate w"^. We 
have 

w'^{y,t) = w^'^ + w^'^ 

where 

w^^\y,t) = 



= -j j G\y + dke3,z + dke3,t~T){diYF + Vp^^''^)iz,T)dzdT 
and 

w''\y,t) = 



= -j j G^{y + dke3,z + dke3,t-T){divF + Vp^^'''^){z,T)dzdT 

In fact, velocity w^'^ obeys the same estimate as (14.11) 

\w^'\y,t)\<c^^. (4.2) 
1/3 + dk 

To prove (14. 2p . we need the following lemma. 
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Lemma 4.1. Let 




Vi{x,t) = I I K^j^{x,y,t - T)Hkm{y,r)dydT 



R3 



2 ^ 2 

wzt/i Kfj^ = Kj^j^ + Kfj^, where K^,,^{x, z, t) is a linear combination of the 
terms 

dVk 

and Kfg^{x, z,t) is a linear combination of the terms 



Gii{x,y,t)— {y,z)dy. 



dxadxji 7 ' ' dys 

Then 



\V{x,t)\<c\\HU-. 

X3 



Proof By ([23]), we find 

2 

\K\x,y,t)\ < ^- '— -^exp ( - ^ 

t2 {\x - y*\^ + t)2 V t 

For the second term K^, we are going to make use of estimate f l2.1ip . So, we 
liave 

t 

R3_ 

t 

R3_ 

t CO ^ 

<c\\H\\^ / f^^ ^ -T exp ( - ^) + )c;z3cir < 



t oo i oo 
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<c||/f||oo^ — J exp(—cu^)du + J — ^— — arctan (^—^^^dr^. 



To estimate the last integral, we do the following 



t t 



/ — ( arctan f — ^^^c/r = / — arctan f — ^dr < — . 



Lemma [4.11 is proved. □ 

Now, upper bound fl4.2l) follows from Lemma 14.11 and the identity 



wf^iy, ^) = y j ^ijmiy + 4e3, z + dkCs, t - T)Fjm{z, T)dzdr 

A IR2x{y3>_dfe} 

with the potential derived from the boundary value problem fl2.7p . in 
which G is replaced with G^. 

If we let w^{y, t) = w^'^ + then the new function is a solution to the 
following initial boundary value problem 



dtw^ - Aw° = div F + Vp 



i(fc) 



in 



w'>{y',-dk,t) = 

for {y',t) G M^x]A,0[, and 

w%y,A) = w{y,A) 

for 1/ G X {t/3 > —dk}- Using exact representation formulae for w^'^ and 
w^'^, we may assume that is bounded by a constant c, which is independent 
of k. On the other hand, we know that function p^^'^), being extend to the 
whole so that p^^'^^y' ,ys,t) = p^'^'^\y' ,ys + 2dk,t) for y^ < —dk, belongs to 
the space Loo{A, 0; BMO{M.^)) and the corresponding norm is bounded by a 
constant independent of k. So, we have 

sup / \p'^'\y,t) - [p'^'%,,,^^it)rdy < c(m)(-A)i. 

A<i<0 J 

-B(x,v^) 
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This means, see Appendix II, that sequence w° is precompact in C{K x 
[A/2,0]), where K is an arbitrary compact of M? . Now, it remains to make 
use of estimates for w^'"^ and w;^'^, pass to the hmit in the equation for 
taking into account that F = w ® w and conclude that, by arbitrariness 
of A, u is Si mild bounded ancient solution to the Navier-Stokes equations 
satisfying |it(0,0)| = 1. □ 

Remark 4.2. In fact we have shown that 

uniformly of on the closure of the set Q{R) for any R > 0. 
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5 Scenario 2 

Here, we are going to prove the following statement. 
Proposition 5.1. 

uniformly of on the closure of the set Q^{R) = B^{R)x] — R'^,0[ with 
B^{R) = {x e B{R) : X3 > 0} for any R > 0. The limit function u is 
equal to zero at — and is not trivial in the sense 

\u{0,a,0) \ = 1 

with a > 0. 

Proof As in the previous section, let us split w into two parts 

w — + 

where 

dtw^ - Aw^ + Vr^ = 0, div = 

in Qa^^Ix]A,01 

w\y\0,t)^0 

for {y',t) e R'^x]A,0[, and 

w\y,A) = w{y,A) 

for y e M^. 

The second part is a solution to the following problem 

dtw'^ - Aw^ + Vr^ = div F + Wp^^''^ , div = 

in Qa, 

w\y',0,t) = 

for {y',t) e M2x]^^o[, and 

w\y,A) = 

for y e Mi^. Here, F = w ® w and p^^'^'' is defined by the following Neumann 
boundary value problem 

Api('=)(y,i) = -divdivF(|/,t), {y,t) e Qa, 
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pf\y',0,t) = {y',t)eR'x]A,0[. 

Let us first discuss precompactness of w"^ in C{B^{R) x [A/2, 0]) for any 
positive R. Indeed, according to (12. Sp . we have 



< 



exp 



t 



and, using estimates fl2.9p and the definition of the kernel K, it is not difficult 
to show that 



j \Kismix,z,t)\dzdT < 



Next, first, assuming A < ti < t2 < and x^,x^ G M.^, 

j j {K{x\ y, h-r)- K{x\ y, h - r))F(|/, r)dydr 



< 



+ 



t2 



+ 




K{x^, y, t2 - T)F{y, T)dydT 



h + h. 



For the second term, we have 



h < c\\F\\ooVt2 - h- 
Next, in the first term, we shall do the change of variables: 



ti-A 



/i<||i^||oo J j \K{x\y,T)-K{x\yM-ti + r)\dydT< 



< \\F\ 




\K{x^, y, t)-K{x^, y, t2-ti + T)\dydT = \\F\\ooJ{x^, x^; ti, ^2 



Using the above estimates of the Green function, it is not difficult to show 
that given e > there exists K,{e, A, R) such that if |a;^ — a;^| < /t, < t2~^i < 
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K, x^x^ G B^{R), and A < ti < t2 < 0, then J(x\ x^; ti, ^2) < So, the 
required precompactness of w"^ follows from Arzela-Ascoli theorem and from 
the following inequality 

\w^{x\ti) -w'^{x'^,t2)\ < \\F\\^J{x\x'^;ti,t2) + c\\F\\ooVh-ti. 

Precompactness of on the same sets is based on similar arguments and 
the following fact 



j \G^ix,y,t)\dy<c 

So, as usual, applying the diagonal Cantor procedure, we select a subse- 
quence still denoted by u^'^^ that converges uniformly on sets Q+(n) for any 
natural n. On the other hand, u^''\0, dk, 0) — i- u{0, a, 0) so that \u{0, a, 0)| = 
1. The latter actually implies that a > 0. □ 

The last arguments in the proof of Proposition 15.11 allow us to prove 
Proposition 11.51 

Proof of Proposition ll.Si If we assume that T is a blow up time, 
then since that a < A we have scenario 2 if one blows up our solution v 
and gets a function u of Proposition 15. 1[ It is not so difficult to see that 
modulus of continuity of u in Q+(2) depends only on the integral modulus of 
continuity of the above Green's functions. So, there exists Sq > such that 
|m(0, 0, 0) - m(0, a, 0)1 < 1/2 provided a < min(l, ^o) but in fact |u(0, 0, 0) - 
M(0,a, 0)1 = |u(0,a, 0)1 = 1. So, it remains to take e = 60. □ 

Proposition 5.2. The limit function u is a mild bounded ancient solution. 

Proof We have already proven that 

u^''^ u (5.1) 

in C(Q4R)). 

We define p^^'^^ as a solution to the following Neumann boundary value 
problem 

Ap^^''\x,t) = -divdiv (^M^^) ® M(^))(x,t) 

for (x, t) G Qt and 

pf\x',0,t) = 
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for (x', t) G M^x] — oo, 0[. We extend p^^'^^ to the whole space in the even 
way with respect to X3 and u^''^ is supposed to be extended by zero. Then 
the function H'^''^ := m^^) ® u^'^) + Ip^^^) G 0; 5M0(R3)) so that 

sup{\\H^^^\\l^(^bmo) + h^'^^Woo) = d<oo. 

k 

We know 

t 

w(^)= I Gix,y,t-A)u^'\y,A)dy+ [ [ K{x,y,t ~ T)F('\y,T)dydT, 



A R3_ 



where F^''^ = u^^'^ ® u^''\ All the norms, bounded by ||u*''^-'||oo, H-^'-'^^Hoo, or 
W H ^''^ \\l 00 (B MO) only, remain to be bounded for the limit functions. 

We may transform the above formula by integration by parts in the second 
term on the right hand side 

wf^(x,t)= [ G,,ix,y,t-A)uf\y,A)dy+ 



dG 




'x, y, t - r){Hf{y, r) - [H^^\^^^,. ,o),a){r))dydr 



dyi 

A ir3 

where B{{x\ 0), a) is a ball of radius a = (xl + t — r)^ centered at the point 
(x',0). 

Then we split u^''^ = u^^^ + u^^^ according to the decomposition of the 
Green function: G = G^ + G'^, see (12. 4p . And let u^^^ = m^^^^ + m^^^ so that 

^mi^^ ^') = f [ ^' ^ ~ ^)(Hj'l\y^ ^) - [Hj'^^]B{ix',0),a)i^))dydT. 



A 



From (ESI), it follows 



dx'^dy'f 



hi 



——{x,y,t-A) <c(a,7)(t-A) ^{t-A + xD 
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Now, we have 



y\x' -yf + yl + xl + t-A)~l. (5.2) 



|Vl"lnJ,5(x, t)\ < c(a)(t -A + x^)-^ x 
X fiW-yf + yl + xl + t- A)-l exp ( - 



and 



A 

X jilx'-yf + yl + xl + t- r)-l exp ( - ^))\H^'\y, r)- 

-[iy('=)]B((^/,0),a)(T)Ml/rfr. 

So, we find 

|Vl"lu(^)(x,t)| < c(a)rf(x2 + t-A)~^ < c{a)d{xl - A)-^-^ (5.3) 
for all (x,t) G ]R^x]A/2,0[ and, by Lemma IGTT] with /3 = ^ and = 



+ t - r. 



|Vl"lM?^^(x,t)| <cia)d {t-TY-^{xl + t-T)--dT < 



< cia)d I < c[a)d I (5.4) 

J {xi + g2) 2 J (xi + g2) 2 

for all (a;,t) G M^x]v4,0[. 

Next, we observe that m^^^ is a solution to the following initial boundary 
value problem: 

9,«;,)-A«;,) = -divif('^) (5.5) 

for any (x,t) G M^x]A,0[, 

n;,)(x',0,t) = (5.6) 
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for any from R^x]A, 0[, 

u\^){x,A) = u^''\x,A) (5.7) 

for all x G Ml. 



Regarding t), we can say the following: for an appropriate exten- 



sion of it to the whole M'^, we can claim that Vmj^^^ is Loo(A, 0; BMO) and 



VM(^) e Loo(A, 0; Lg,unif(IR^)) for any finite g > 1. 

Since H^^'^ is bounded in Loo (A, 0; BMOiB?)), by estimates O and ^ 
we have (from the energy inequality) 








sup j I (p'^l'Vuj^f.-^l'^dxdt < oo (5.8) 

A/21 



for any test function G C(j"(M'^x] — A/2, A/2[). Then, from the equation 
for the pressure p^^''\ we derive the similar estimate for the gradient of p^^^-* 



sup y j ifi'^\Vp^^''^\'^dxdt < oo (5.9) 

A/2RI 

for any test function (p G C^(M^x] — A/2, A/2[). Using coercive estimates 
for the heat equation, we then get 



sup J y v9^(|9tU(fc)p + |V\(fc)nc/xc/t < oo (5.10) 

A/21R3 

for any test function <p G ^^^(M^x] - A/2, A/2[). 

Next, we may use a parabolic imbedding theorem to show that 



sup J J if'^l'Vu^^i^^fdxdt < oo 

A/2R3 

for the same test function (p. The same type of estimate is valid for the 
gradient of p^^*^^ . Repeating these arguments several times, we can state that 
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Vm^j.) is bounded in domains M"^x]3A/8, 0[. Boundedness of Vp^'-'^^ follows 



from local regularity and the equation Ap^^^^ = —u'i'j'uf} being valid in 



for appropriate extensions of p^^^^ and uf^^u]} from Mi to M^. Summarizing 



all these estimates, taking into account arbitrariness of A, and using a shift 
in time we find 

(5.11) 



sup \Vu^^\x,t)\ < cln(2 



for (x, t) G Qt and 



sup sup \Vp^^''\x,t)\ <C<oo. 

k (x,t)eQl 



(5.12) 



Now, we can go back to evaluation of function u^^^. By fl5.1ip and fl5.12p . 



we do not need integrate by parts in the expression for m^^^ any more and 
thus 



\{k) 




d 



GUx,y,t-r) ^Ff{y,T) + ^p'^'\y,r) dydr. 



9 1 



\-dyi 



Next, according to (15. lip and (I5.12p and by estimates of Green function 



we find the following bound 

t 



|VM(')(x,t)| < C 



dr 



xl + t-T)U {xl + yi + t-T[ 



-ln(2 + l/?/3)c?l/3 



< c 



dr 



exp - ^ 



xl + t-r)-^ J {xl + yl + t-r)-i 



-ln(2 + l/|/3)%3 



< c 



dr 



(xi + t-r)4 



^3 2 ln(2 + l/y3)dy3 + / exp 



t - T 



dy-s 



< c(i + v^)(-A)i 

The latter implies that in fact 

sup sup |VM^''^(x,t)| < C < 00. 



(5.13) 
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Using estimates f l2.8p . fl5.2l) . fl5.12p . and bootstrap arguments, we show 
sup |V'M(^)(x,t)| + |V'+V('^Hx,t)| < 0(1,5) < oo (5.14) 

for any / > and for any 6 > 0. 

As to derivatives in t, by fl2.10p . we find 



t oo 



So, 




A 



■X 



X exp I — 



< c 



t - T 



divif(^')(y,r) 
1 



(t - r)2 (xl + t - r)2 



sup \dtuf^Jx',X3,t)\ <c [ 



1 1 



2 (x| + 13)2 



For the first term, we show in the same way that 

\dtu^^^{x,t)\<c I \dtG\x,y,t-A)\\u^^\y,A)\dy 



(5.15) 



and thus 



\dtU(k){x,t)\ < 



(t - A)^ (xj + t- A)^' 
From equation (15. 5p . (15.150 . and (I5.16p . we can deduce that: 



sup \dtu'^''\x,t)\ <c{6). 

{x,t)&Qt,xs>S 



(5.16) 



(5.17) 



After passing to the hmit as A; — )■ oo, we get that u satisfies the standard 
integral identity with divergence free test functions. Then one can claim that 
there exists such that 



dtu — Au + Vj9^ = — div H 



(5.18) 
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in gl. 

Obviously, p'^{-,t) is a harmonic function in a half space whose gradient 
is bounded in t, in x', and in > 6 for any 6 > 0. Taking the limit in f lS.Sp . 
we show that Vp^ satisfies 

dtu^ - Au^ + Vp^ = 
in Mi^x]A,0[. This, together with O, © and f lSlSjl . ([EISD, implies 

sup |Vp^(a;', xs, t)| — )■ 

as rca ^ CX3 for each — oo < t < 0. Moreover, more detail analysis of the 
above estimates for Green functions allows us to state: 

|Vp2(x,t)| < cln(2 + I/X3) 

for all {x,t) G Q^. Proposition 15.21 is proved. □ 

Proof of Theorem II .41 easily follows from the above arguments. 
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6 Appendix I 



Lemma 6.1. Assume that numbers m and satisfy the condition 

< ao < (6.1) 
Let a and P be positive numbers and let us define 

T( R\ f \f{x) -[f]B(a)\ ( Xl\ 

where [f]B{R) is a mean value of f over the ball B{R) of radius R centered 
at the origin. Let further n — fUi — j^^, and rii — J^^^ . Then 

3qq 

|/(a,/3)| < c(m,Q;o)||/||BMO(R3)a" - x 

oo 

X (v^ / exp(-«^) ^Y. (6.2) 

Proof Without loss of generality, we may assume that [f\B{a) — 0. We 
also let 

1 / \ 

^(a;) := — ^ — ;;^^^p(~^)' ^ ■= \\f\\BMO{M?)- 

Then by Holder inequality we have 

\I{a,l3)\< J K^^^{x)\f{x)\dx < 

<(^J K^{x)\f{x)rdxY j K^^{x)dxy. 

R3 ]R3 

For the second factor on the right hand side of the latter inequality, we find 

oo oo 

jK^.ix)dx^ Je^{-^)dx,2.l 
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oo 

Iff xhi . , 
2^-3;i 7 / exp ( - 7^ ]dx^ 



2ni 



-I J V pmJ u2.^2\2^^-^ 







Here, we have used condition (16 .ip that ensures the inequahty ^ — 1 > 0. 
After changing variable, we have 

^ oo ^ 

K~{x)dxY <c{m,a)(^ / exp s;^)"- 

^ ^ J {u'^Pni/n + a?)^~^^ 

Next, 

/m I m 

K~i{x)\f{x)rdx = / i^^(x)i/(x)rrfx+ 

R3 B(a) 
oo „ 

+ / ir^(a;)|/(x)rda; = /i + /2. 

''=°B{a2'=+i)\_B(a2'=) 

Obviously, 

For J2, we are going to use the fact that m/mi = 1 + ao and two inequalities 
and 

|[/]s(a2^+i)| <c(A; + l)A 

The latter inequality can be found for example in [15], p. 141. Then, we 
have 

B(a2'=+i) 

00 ^ 

fc=0 ^ 

This completes the proof of Lemma I6.1[ □ 



27 



7 Appendix II 



Consider the following problem 



dfU — Au — — div g 



in Q{a) = Q{a) + (0, a^). Assuming that u and g are sufficiently smooth in 
Q{a), we wish to estimate the modulus of continuity of u in the closure of 
the set B(a/2)x]3(a/2)2,a2[ in terms of ||ii||oo,Q(a) and |b - [5']s(a) IU^,oc(S(a)) 
for sufficiently large m > 1. 

Fist we take a positive cut-off function ip e C^(5(a) x]0, 2a^[) and ex- 
tended it by zero. We let w — (pu. This function solves the following Cauchy 
problem 

dtw — Aw = —div (iff) + fV(p + w{dt(p + A(p) — 2div (mV<^) 

in ]R^x]0,a^[, where f{x,t) = g{x,t) - [g]B{a){t), and w{x,0) = for x e R^. 
We split w into two parts w — + w"^, where 

dtw^ - Aw^ = -div (ipf) - 2div (uV</?) 

in M^xjO, with w^{-, 0) = in and 

dtw'^ - Aw'^ = fV(fi + w{dt(fi + A(p) 

in M^x]0, with 0) = in M^. The second part can be estimated with 
the help of coercive estimates 



For sufficiently large m, they give an estimate for a Holder norm of u and 
thus we have a required estimate for the modulus of continuity of w"^. 
For the first part, let us make use of the solution formula 



w^{x,t)^- / / T {x - y, t-r) [div {if f){y,T) + 2diy{uV(p)){y,T))dydT 




W^i(M3x]0,a2[) 



< C{a,m, \\u\\oo,Q{a}, ||/||L^,oo(Q(a)))- 




R3 




R3 
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We first assume that m > 2012 and observe that the following simple 
estimate is true: 



J \K{x-y,t-T)\dy<-j^, 



where K{x—y, t) = Vx^{x—y, t). Now, we wish to show that a given positive 
s there exists a positive 6{m,a) such that if {x^,ti), {x^,t2) G B{a)x]0,a^[ 
with ^2 > ^1 and \x^ — x^\ + t2 — h < S, then 








\K{x^ -y,t2-ti + ^) - K{x^ - y, i3)\dy) ^ —dd < e. 



Assume that this is false. Then there exists £o > and sequences 

(a;^'",ii,„),(a;2'-,t2,n) e5(a)x]0,a2[ 
with t2,n > ti,n and — + t2,n — but 

M3 

Indeed, for any positive 

by the Lebesgue theorem. And then 

/" ( /" \K{x'''--y,t2,n-h,n + i^)-K{x'^^-y,d)\dy]'^^d^^0, 



rs 
since 



— 1 1 



•dm '{}2m' 



c 

< ^ — 



1?2 + 2rn, 
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Since the function on the right hand side is integrable under our assumption 

on m, we arrive at a contradiction again by the Lebesgue theorem on the 
dominated convergence. The rest goes as follows: 

\w\x'^,t2)-w\x\ti)\ < 

t-z 

J J K{x2-y,t2-T)-(^{ipf){y,T) + 2{uVip)){y,T)yydi 



< 



+ 



tl R3 



+ 



J J {K{x2-y,t2-T)-K{x,,h-T))-(^{ipf){y,T)+2{uVip)){y,T)yydT 



R3 



Then we apply Holder inequality and have 

\w\x'^,t2) -w\x\ti)\ < 



< 




K{x2 -y,t2- t) 



dy) X 



x( j \K{x2-y,t2-r)\{\fr + \ur)dydT) 



+ 



B{a) 



+ 




K{x2 -y,t2-T) - K{xi,ti - t] 



dy)"' X 



R3 



x( J |i^(a;2-2/,t2-r)-i^(a;i,ti-r)|(|/|"+|Mr)rf2/rfr) 



< 



B{a) 



t2 

<C{a,m,\\u\\oo,Q(a)A\f\\LmMQia)))( / t: 7TT^d^+ 

{t2 — T)2^2m 



+ [ {f \K{x^-yM-h + ^)-K{x^-y,myY^d{i) 



- C(a,m, ll-ulloo.Qia), ll/llw(Q(a))) (^2 -ii)^ 2m (i+ 
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R3 

This, together with the above statements, gives a required estimate for the 
modulus of continuity in the closure of the -B(a/2) x]3(a/2)^, o?[. 
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